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, We investigate the threshold behavior of transmission resonances and quasibound states in the 

■ multichannel scattering problems of a one dimensional (ID) time-dependent impurity potential, 

^Sj ' and the related problem of a single impurity in a quasi ID wire. It was claimed before in the 

literature that a quasibound state disappears when a transmission zero collides with the subband 
O ' boundary. However, the transmission line shape, the Friedel sum rule, and the delay time show 

\ that the quasibound states still survive and affect the physical quantities. We discuss the relation 

between threshold behavior of transmission resonances, and quasibound states and their boundary 
conditions in the general context of multichannel scatterings. 
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I. INTRODUCTION 



Electronic transport properties in quasi one-dimensional (ID) wires have gained much attention, not only due to 
d . scientific interest but also due to practical applications [1]. When an electron is restricted to a wire, confinement 
subbands are formed due to quantization of the transverse momentum. In clean straight wires without impurities the 
I . modes of the subbands decouple, which results in a perfectly quantized conductance. However, the modes (including 
evanescent ones) mix when impurities are introduced into the wire, and the transport properties show a rich and 
Q ■ non-trivial behavior. Detailed aspects of the mixing have been extensively investigated [2-7]. The impurities put into 
^ the wire usually give rise to bound states, which can be separated into two categories: true bound states with real 
energy and quasibound states with complex energy. The imaginary part determines the decay rate of the quasibound 
state. 

^ . Even in simple ID scattering problems quasibound states are quite common, for example resonant transmission 
" states through a double barrier [8,9] and Ramsauer-like resonances [10] in potential wells. According to the Breit- 
Wigner theory [11], for each quasibound state the transmission amplitude possesses a pole in the complex energy 
plane. In multichannel scattering problems another kind of quasibound states can appear; the Fano type resonance 
originally proposed to explain autoionization in atomic physics [12]. When a discrete energy level interferes with a 
' continuum of states, the excitation spectra become asymmetric. Each subband mode in a quasi-lD wire gives rise 
. to a zero-pole pair of transmission in the complex energy plane, where the zero occurs on the real energy axis. As a 
consequence, there exist transmission zeros, and the zero-pole pairs lead to asymmetric resonance features [13]. 
Transport through a time-dependent (mostly periodically oscillating) potential is also a subject of increasing im- 
S ' portance, with a growing number of applications [14-21]. One of the important features here is that an oscillating 
. potential can transfer an incoming electron of energy E with finite probability to 'Floquet' sidebands at i5 ± nliio, 
(~| ■ where n is an integer and w is the angular frequency of the oscillation. Even though this corresponds to an inelastic 
O ■ scattering, the process is coherent and non-dissipative. Hence, it can still be regarded as a multichannel scattering 
^ problem like in the quasi-lD case, with the sidebands being analogous to the subbands. However, from the physical 
^ point of view, the Fano type quasibound states of the periodically oscillating impurity are somewhat different from 
those of quasi- ID system [22]. If a bound state exists due to an attractive static potential added to the oscillating 
, potential, electrons in the incident channel can emit photons and drop to the bound state, and similarly electrons in 
d ' the bound state can absorb photons and jump to the incident channel. A transmission resonance takes place when 
the energy difference between the incident channel and the bound state is equal to the energy of one or more photons 
of the oscillating potential. Even without bound states transmission resonances like these can occur. 

The multichannel subbands in the quasi-lD case and the Floquet sidebands in the oscillating potential case divide 
the complex energy plane into sub-domains with appropriate boundary conditions for the quasibound states. Upon 
varying an external parameter such as the strength of the impurity potential, the quasibound states change their 
location in complex energy plane. In the literature [20] it was claimed for the oscillating impurity that the quasibound 
state vanishes when the transmission zero collides with the sub-domain boundary and disappears. However, as we 
will demonstrate, the transmission line shape, the Friedel sum rule, and the delay time show that the quasibound 
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states do not disappear abruptly. This clarifies the relation between threshold behavior of transmission resonances, 
and quasibound states and their boundary conditions in the general context of multichannel scatterings. 

In Sec. II, we briefly summarize our findings about the transmission resonances and quasibound states. In Sec. Ill, 
we investigate the characteristics of quasibound states in the scattering problem with a ID oscillating delta-function 
impurity, and discuss Friedel sum rule and Wigner delay time. In Sec. IV, we study our second example, the scattering 
from a static delta-function impurity in a quasi-lD wire. Finally, we conclude our paper in Sec. V. 



II. BOUNDARY CONDITIONS OF QUASIBOUND STATES IN MULTICHANNEL SCATTERING 

In this section we summarize the main finding of this paper, which will be further illustrated in Sec. Ill and IV. 
The transmission coefficient tmn{E) relates the incoming wave V'JT to the scattered wave by 

rr:i{E)=tmn{E)r^. (1) 

The resonance energies are found by locating the poles of tmn{E), namely 

^ =0. (2) 



tmn{ER — lEl) 

The poles of tmn{E) occur at complex energies E — Eji — iEj. At these complex energies, a scattered wave can be 
produced even if no incident wave is present. When we choose the sign convention of time as exp{—iEt/?i), Ej must 
be positive for the system to be stable. In the Breit- Wigner case [11] the transmission near the quasibound state is 
given by t{E) ^ 1/[E — {En — iEj)], whore the pole in complex energy plane completely determines the transmission 
line shape. 

The poles can also be found by solving the wave equation with the appropriate boundary conditions of only outgoing 
waves. This can be described by the following linear equation 

^Lt \ _ 



Here M is the matrix of wave-matching conditions, and V'iut ^^^d "^otit the left and the right outgoing wave vectors, 
respectively. The matrix M contains the wave number fc as a function of energy E. For real E the wave number k 
is given by fc = ^/2jlE/fi, where /i is the mass. For complex energies we have to choose the sign of the square root 
carefully. Considering the general relation [23] 



provided that /3 ^ 0, we obtain 



(4) 



VEh - iEi = ±{kR - iki), (5) 

where by definition Ej, kn and hi are positive. The real energy E^ > for propagating states, and Efi < for 
evanescent states. The correct choice of the sign in Eq. (5) depends on the analytical continuation based on the 
physical constraints on the real energy axis. For Ej = the following relation should be satisfied: 

where both ka and kj are still positive. The choice of sign for propagating states corresponds to the choice of the 
negative sign for kj (i.e., the positive sign for ka), and follows from the propagating direction of the plane wave far 
from the scatterer. For evanescent states the positive sign for ki must be chosen in order to avoid divergence of the 
solution in space. 

It is worth noting that the '-h' sign corresponds to the decaying wave exp{ikRX — kjx) exp{—iEiit/h — Ejt/h) while 
the ' sign corresponds to the diverging wave cxp{ikRX + kjx) exp{—iEjit /h — Eit/h), with a; > 0. The physical 
reason of the exponential divergence of the channels with positive real energy is a retardation effect; the wave at 
X ^ 1 has propagated away from the impurity where it originated at a time At w x/v in the past, where v is the 
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velocity of the wave. However, at that earher time the anipHtude of the wave at the impurity was larger by a factor 
exp{EiAt/h). This corresponds to exp(A;/a;) since {Er — iEi)/h ~ v(kR — iki)/2 [24]. 
For multichannel scattering the wave numbers are given by 

kn = ±^/Er -En- iEi, (7) 

where En represents the quantized energy of a transverse mode in the quasi-lD case (n = 1, 2, 3, • • •) and is equal to 
nUw in the ID oscillating potential case (n = • • • , —1, 0, 1, • • •). We order the E^s by magnitude and denote by n* 
the special value of n which satisfies < Er < En'+i- 

In order to satisfy the physical conditions for the signs of A;„'s mentioned above, for n < n* the '— ' sign should 
be chosen for Imfc„, and for n > n* the '+' sign, which corresponds to the selection of a certain Riemann sheet for 
evaluating the complex square root function. This choice of the signs will be called proper Riemann sheet since this 
set of signs corresponds to the correct analytical continuation from the real energy axis. In order to represent the 
Riemann sheet through the signs of Imfc like in our convention, the square root functions should have their branch 
cuts on the real axis starting from the branch points and extending to the positive infinity. In general the notion 
Riemann sheet depends on the choice of a branch cut. However, this choice of the branch cut cannot affect physical 
results, because it is only a way to label the projection of the Riemann sheets onto the complex plane. 

We introduce a labeling for the Riemann sheets by the vector of signs of the square root function. In this notation 
the proper Riemann sheet is represented by (•••,—,—, G, +,+,•• •), which implies '— ' for Imfc„*_2, '— ' for ImA;„*_i, 
'— ' for Imfc„. , '+' for ImA;„*+i, and '+' for Imfc„*+2, and 9 has been used to indicate the sign of n*. We also call the 
regions of complex energies separated by the quantized energies En's "domains", e.g. Z)„ represents the region with 
En < Er < En+i as shown in Fig. 1(a). 

As a parameter of the scattering potential such as its strength is varied, the locations of the poles change continuously 
and sometimes can collide with the boundary of the initial domain. We can trace the pole into the neighboring 
domain through the boundary, for example from D„. to Dn*-i in Fig. 1(a). In fact the trajectory of the pole 
across the boundary is continuous and smooth once we use the same Riemann sheet as the used one in D„. , i.e. 
(•••,—,—, 0, +,+,•• •). However, this does not correspond to the proper Riemann sheet (•••,—,—,©,+,+,•••) in 
this new domain Dn*-i since the inequality En* < Er is violated. Thus, the solution obtained in Dn*-i using the 
proper Riemann sheet of Z)„. does not have any physical meaning in the domain Dn*-i (more precisely on the real 
axis o/_D„*_i), but still can influence physics on the real axis of -D„*, especially when the resonance is broad. Let 
us emphasize that only the branch points. En's, play an important role in the threshold behavior of transmission 
resonances, but not the placement of the cuts. Physics is in the choice of signs on the real energy axis, Eq. (5). This 
situation is summarized schematically in Fig. 1(b), which shows that only small part of the proper Riemann sheet of 
Dn*, i.e. En* < Er < En*+i, is attached to the real energy axis. Since in either Dn*-i or Dn*+i the quasibound 
state can survive on the proper Riemann sheet of Z)„* , it can still affect the physics in the region En* < Er < En*+i, 
but has no effect on other parts of real axis. In the following sections we discuss two examples. 

III. ID DELTA-FUNCTION IMPURITY WITH TIME DEPENDENCE 

As the first example we investigate a ID scattering problem with a delta-function impurity oscillating at frequency 
Lu. This problem was studied before by Bagwell and Lake [4], who calculated electron transmission and discussed 
transmission resonances. Recently Martinez and Reichl [20] have investigated this problem and found the existence of 
nonresonant bands in the transmission amplitude as a function of the strength of the potential and driving frequency. 
They observed a periodic behavior of the conductance as a function of the scattering strength and of the oscillation 
frequency of the scatterer. 

A. Scattering matrix formulation 

The system is described by the Hamiltonian 

H{x,t) = -^-^ + [V, + VdCos{iot)]S{x), (8) 

where the fi is the mass of an incident particle, and Vs and Vd represent the strength of the static and the oscillating 
potential, respectively. Using the Floquet formalism [25] the solution of this Hamiltonian can be expressed as 
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x)e 



-inujt 



(9) 



where Epi is the Floquet energy whic;h take continuous vahies in the interval < Epi < 'huj. 
Since the potential is zero everywhere except at a: = 0, ipn{x) is given by the following form 



-iknX 



x<0 
X > 0, 



where fc„ = \/2ii{Epi + nTiuj)/%. The wave function ^EFi{x,t) is continuous at a; = 0, 
and the derivative jumps by 



dx 



i;=0+ 



dx 



2m 

= -^[ys + VdCos{u;t)]^EF,{0,t)- 



Using Eq. (9) this leads to the condition 

ikn{Cn D 



An + Bn) = ls{An + B„) + 7d(A„+i + A„_i + B„+i + B„_i) 

= 7,(C„ + AO + 7<i(C„+i + C„_i + £»„+! + 



(10) 



(11) 



(12) 



(13) 
(14) 



where 7^ = 2ijlVs/Tt? and 7^ = iJ,Vd/h^. After some algebra we have the following equation from Eqs. (11), (13), and 

(14) 



-(/ + r)-ir (/ + r)-i 
(7 + r)-i -(7 + r)-ir 



(15) 



where 



r = 



-.000 

jd/ik-1 Is/ik-i Id/ik-i 
Id/iko Is/iko Id/iko 



V 



^d/iki ^s/iki Jd/iki 

■•• 



(16) 



and / is an infinite-dimensional square identity matrix. Eq. (15) can also be expressed in the form |out) = M |in), 
where M connects the input coefficients to the output coefficients including the associated evanescent Floquet side- 
bands. In order to construct the scattering matrix we multiply an identity to both sides, K~^K |out) = MK~^K |in), 
where Knm = \/kn^nm- Then we have Jout = MJin, where J represents the amplitude of probability flux and 
M = KMK"^. It should be mentioned that M is not unitary due to the evanescent modes included. 

If we keep only the propagating modes, we obtain the unitary scattering matrix S [21,22], which can be expressed 
in the following form 



(17) 
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where r„m and tnm are the reflection and the transmission amplitudes, respectively, for modes incident from the left; 

milar quantities for modes incident from the r: 
coefficient for the propagating mode entering in the mth channel 



r'nm and t'^^ are similar quantities for modes incident from the right. From S we can obtain the total transmission 
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oc 
n=0 



(18) 



Figure 2 shows the transmission cocfScicnt T as a function of s for various values of a dimensionless potential 
strength a^, where e = E/Tuv [E \s the kinetic energy of the incident particle) and aa = /iV^/Sh^u;. We exploit the 
relation T{E) = TmiEpi) where E = Epi + mTvjj. For small o^'s the transmission coefficient shows some signatures 
of Fano type transmission resonance structures [Fig. 2(a) and (b)] whereas the resonance structure and transmission 
zero disappear as we increase [Fig. 2(c)]. Although there exists no transmission zero in Fig. 2(c), the resonance- 
like structure (the peak around e = 0.8) still remains. For much larger the transmission zero reappears from 
e = 1. Let us note that the transmission coefficient has no symmetry of energy translation by mTiuj since in general 
{Epi) ^ Tm'iEpi) when m^m'. 



B. Quasibound states 



In order to obtain quasibound states we apply the same procedure used in the previous subsection, but with the 



different boundary conditions, A = Q and D 
little algebra we obtain the linear equation 



0. Since i?„ = C„ because of the reflection symmetry at a; = 0, after 
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V-) 



(19) 



Here = iJ.V^/2n^uj. We solve Eq. (19) for given ag and by using singular value decomposition [26] and find the 
complex energies of the quasibound states. 

Due to the structure of the matrix in Eq. (19), if e' is a solution of Eq. (19), then e' + n is also a solution, where n 
is an integer. We limit ourselves to the energy range of < < 1, which means the poles are meaningful only in the 
domain Do- This constraint determines the signs of the complex square root functions in the diagonal of the matrix 
in Eq. (19) as described in Sec. II. Using the same notation as the one used in Sec. II the signs can be expressed by 

(•••, +,+,e, -,-,•••)• 

Figure 3(a) shows a trajectory of the quasiboimd state as is varied with a, = 0. The quasiboimd state of the 
first branch from the top originates from (1,0) at ad = 0, moves towards the domain -D-i, and crosses the domain 
boundary er = ^ sX ad — 0.935. The trajectory of the pole helps to understand the transmission curves in Fig. 2. The 
resonance structure of the transmission follows the pole, as is shown in Fig. 2 and the inset of Fig. 2(a). Since in the 
system considered here the transmission zero always takes place on the left of the pole in the real energy direction, 
it collides with the domain boundary er = 1 prior to the pole and disappears first at — 0.782. In Ref. [20] it 
is concluded that the quasibound state vanishes when the transmission zero disappears. Even after the pole moves 
out of the domain Do, however, it still affects the transmission line shape, as shown in Fig. 2(c). There is no abrupt 
change of physical quantities, e.g. the transmissions, the Friedel phases, and the Wigner delay time (these two will be 
shown below). The transmission zero disappears when it collides with the point (0, 0). We can trace the transmission 
zero even into the non-physical region er < Q hy calculating S matrix on the proper Riemann sheet of Dqj although 
the unitarity of the S matrix is violated. 

A quasibound state of a second branch appears from the real axis at ~ 1.032 and moves to the left. It passes 
through the domain boundaries er = 1 &i ad — 1-31, while the transmission zero appears in advance at ad — 1.165 
as shown in Fig. 2(d). The pole finally leaves the domain Dq by passing = at ad — 1.935. The second branch 
of the pole exactly corresponds to the case of the schematic shown in Fig. 1(b). In fact one can find more branches 
of quasibound states in a regular fashion [20]. It is noted, however, that the branches qualitatively differ in their 
solution vectors, i.e., the null space of the matrix of Eq. (19) at the energy of quasibound state as shown in the insets 
of Fig. 3(a). 

Figure 3(b) shows the trajectory of the pole in the case of a delta- function impurity with an attractive static 
potential added to the oscillating potential. A similar behavior as in Fig. 3(a) is observed. The shifted starting point 
from £ = 1 at ad = is ascribed to the existence of the true bound state of an attractive delta function potential with 
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binding energy Eb = — /iV//27i^. The introduction of the static potential makes no effect on the signs of the complex 
square root functions, i.e., the Riemann sheet, which is clear in Eq. (19). 

An interesting difference takes place in the trajectory of the quasibound states in the oscillating delta function 
impurity when the attractive static potential is present. For small the matrix in Eq. (19) can be approximated by 
a 3 X 3 matrix involving only k-i, ko and fci, which corresponds to single-photon processes (absorption and emission). 
Using the results of Appendix B of Ref. [4] it can be seen that this 3x3 matrix is not enough for the case without a 
static potential to describe the shift along the real energy axis of the quasiboimd states as is increased. At least 
a 5 X 5 matrix is needed to include the shift, which implies two-photon process. The mathematical reason is the fact 
that the slope dei/deala^^o = oo in Fig. 3(a). On the other hand, in the case with an attractive potential 3x3 
matrix is enough to obtain a shift of the pole along the real energy axis, as shown in Fig. 3(b) because dej/den\a^=o 
has a finite value. With an attractive potential, a bound state already exists, and the oscillating potential leads to 
an ac-Stark shift [27] of that energy level. Without a static potential, however, one photon is needed to generate a 
bound state, and additional photons are needed to produce the ac-Stark-like shift. 

C. Friedel sum rule 

The number Nb of bound states below an energy Ef (say, the Fermi energy) is related to the Friedel phase 
6p = lnDet{S)/2i of the scattering matrix S though the Friedel sum rule [28] (see also [29]) 

Nb = -OpiEf). (20) 

TT 

This implies that each bound state contributes tt to the Friedel phase. Let us note that the scattering matrix in 
Eq. (17) is given as a function of a Floquet energy epi (= Epi/TiLo) and the same holds for the Friedel phase, i.e. 
it is only relevant for < e < 1. Figures 4(a) and (b) show Friedel phases as a function of epi for = 0.6 (the 
pole is located in Dq) and = 1-0 (the pole is in f-i), respectively, with ag = 0. No abrupt change of the Friedel 
phase is observed at the threshold aa — 0.782 of the transmission zero or when the pole crosses the domain boundary, 
Od ~ 0.935. 

It is noted in Fig. 4(a) that the Friedel phase does not exhibit a clear tt shift but looks even similar to Fig. 4(b). 

We expect the Friedel phase to shift by tt near quasibound states. The absence of a clear n shift in Fig. 4(a) is 
mainly ascribed to the large overlap between the resonance and the neighboring domain since the imaginary part of 
the complex energy of a pole ej is quite large, and so is the width of the resonance. Even though sj is quite small for 
very small a^, the overlap is still large enough to violate the tt phase shift by Friedel sum rule as shown in Fig. 4(c). 
The quasibound state is located quite close to the boundary = 1 because of dei/deii\a^=o = oo. If we consider the 
oscillating potential with an additional attractive static potential, however, for very small aa a gap exists between the 
bound state and the domain boundary er = I, which can make the overlap negligibly small, so that the phase shift 
TT is clearly seen in Fig. 4(d). When the Friedel sum rule is applied to the case of a multichannel scattering problem, 
the overlap between the resonance and neighboring domains should be carefully considered. 

D. Wigner delay time 

The time scales associated with the quantum scattering have been a quite controversial issue (see e.g. [30]) because 
time is not a Hermitian operator in Hilbert space but a parameter. Nevertheless one can define the global time delay 
of wave-packets by using a Hermitian time delay operator (see e.g. [31]). This time delay operator turns out to be 
related to the so-called Wigner-Smith time delay matrix [32], whose trace averaged over different channels gives a 
measure of the time delay of the scattered wave caused by the scattering potential (see also [33,34]). The Wigner 
delay time shows a peak at resonance energies, where the maximal value of the Wigner delay time approximately 
corresponds to the decay time h/ej. To obtain the Wigner delay time we use the eigenvalues of the scattering matrix 
S. Due to the unitarity of S, all the eigenvalues lie on the unit circle and can be written in the form exp{i6a)- The 
Wigner delay time is defined by 

rw = nY^^\{kM\\ (21) 

a 

where the eigenstate corresponding to the eigenvalue 9a and an input propagating state (or channel) with momentum 
kn are denoted by \9a) and |A;„), respectively [35]. It is worth noting that the Wigner delay time is a function of the 
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energy of the incident particle. The insets of Fig 5 shows the Wigner delay time as a function of energy of an incident 
particle for two a^'s with = 0. The peaks of the Wigner delay time are located almost exactly at the real part 
of the complex energy of the pole in Fig. 5 (We omitted the data close to the domain boundaries where the Wigner 
delay time diverges), imder the condition that the pole is located in Dq. No abrupt change of the Wigner delay time 
is again observed at the threshold of the transmission zero (a^ = 0.782). However, the peak disappears close to the 
domain boundary Z)_i since the pole enters the next domain. 



IV. DELTA-FUNCTION IMPURITY IN QUASI-ID WIRE 



As our second illustrative example we study the simple quasi-lD scattering problem with an attractive static delta- 
function impurity of strength |7|. In Ref. [3] it was found that dips in the transmission disappear as |7| is increased. 
Let us again discuss the relation to transmission poles. 



A. Scattering matrix formulation 

The Hamiltonian is given by 

H = j^+l5{x)S{y-yo) + Ve{y), (22) 

where 7 < 0, and Vc represents a confinement potential which restricts the movement of the particle to the range 
< y < D in transversal direction. Using the complete set of the transversal modes xiv) = \/2/£> sm{mTy /D), the 
solution of this Hamiltonian can be expressed as 

00 

i;E{x,y) = J2('n{x)x{y)- (23) 

n=l 

Plugging this expansion into the Schrodinger equation and exploiting orthogonality of the transversal modes, we 
obtain the following equation 

+ klCmix) = ^ MmnCn{x)d{x), (24) 



where Mm„ = {Afij/ Dh^) sm{imryo/ D) s'm{nnyo/ D) and kn = \J2^E/Ti^ — (mr/D)^. Using both the plane waves of 
Eq. (10) for c„(x) and the constraints in Eqs. (11) and (12), we obtain the following conditions 

A„ + S„=C„ + D„, (25) 
ikniCn - Dn) - ikn{An - B„) = ^M„„(A„ + Bm)- (26) 



These can be rewritten as 



B \ _ ( {2Q- M)-'^M 2{2Q - M)-^Q \ ( A 
C )~\ 2{2Q - M)-^Q (2Q - M)-^M )\d 



(27) 



where Qmn = ik.a8.„in- Following the same procedure as that used in Eq. (17) we can obtain the ^-matrix for this 
problem and the total transmission coefficient 

mn(prop.) 

where the sum extends over all propagating modes. The resulting transmission coefiicients is shown in Fig. 6. We 
use the parameters D = 30 nm for the width of the wire, the mass n = 0.067me as the efi'ective mass of an electron 
in GaAs, yo = 5D/12 for the transversal position of the delta- function scatterer, and energy E normalized by the 
El = 6.24 meV. The transmission coeSicient depends on the number of the propagating channels, which will be 
discussed in the next subsection. 
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B. QuEtsibound states 



Quasibound states can be obtained from the following equation, derived from Green's frmction approach [36] , 



which has to be solved for energy Eji — iEj [entering via Eq. (7)]. It should be mentioned that the location of the 
quasibound states and the transmission coefficients strongly depend on the number of modes Nch [36]. In fact they 
do not converge although Nch is increased. In order to model a realistic impurity one should take Nch large but finite, 
and fix it for all numerical calculations. N^h = 100 is chosen in this paper. 

Figure 7 shows the trajectory of a quasibound state as 7 is varied. The quasibound state crosses the domain 
boundary £^2 at 7 ~ —7.11 feVcm^, and enters the next domain D2 (the open circles in Fig. 7) on the proper Riemann 
sheet of the domain Di, (Q, +,+,•• •). The situation is similar to those out of Sec. Ill except that the quasibound 
state collides with the real axis at 7 ~ —9.06 feVcm^ and only one branch of quasibound states exists. As I7I increases 
the transmission zero moves far away from the pole, as is shown in the inset of Fig. 7, and finally vanishes when it 
collides with the other domain boundary i^i at 7 = —8.83 feVcm^. 

Apart from the quasibound state, the attractive delta-function impurity also accommodates one true bound state 
below El . Although the quasibound state vanishes at the threshold value of 7, the true bound state exists regardless 
of the value of 7, and its energy monotonically decreases as I7I increases. 



The Friedel phase of the impurity in the quasi-lD wire is shown in Fig 8. The tt phase shift by quasibound state is 
clearly observed for very small I7I, but it is hard to be resolved for larger I7I due to the overlap between the resonance 
and the neighboring domain. Figure 9 shows the peaks of Wigner delay time and the real part of the complex energy 
of the pole, which coincide very well for I7I < 5.5 feVcm^. For larger I7I it is again difficult to distinguish the peaks of 
the Wigner delay time from the background which diverges at the domain boundary E2 (see Fig. 7). As I7I increases 
the pole first moves away from E2, and then moves back toward E2. The peak of the Wigner delay time originating 
from the pole strongly overlaps with the divergent part of the background when it changes the direction. Hence, for 
large I7I the position of the peak no longer coincides with the energy of the pole. Friedel phase and Wigner delay 
time of the quasi- ID case behave similarly to the case of an oscillating impurity. 



We have investigated the characteristics of transmission resonances, Friedel phases, and Wigner delay times in 
multichannel scattering problems focusing on the trajectories of quasibound states in complex energy plane as a 
scattering strength is changed. Subband thresholds divide the complex energy plane into sub-domains characterized 
by appropriate boundary conditions for the quasibound states. This set of boundary condition is related to the choice 
of the signs of the complex wave number fc„ in each subband. A transmission resonance (or, a transmission zero) can 
vanish when it collides with the subband boundary, which, however, does not mean that the pole becomes completely 
irrelevant, as is seen in the transmission line shape, the Friedel phase, and the delay time. Even if a quasibound state 
enters the neighboring domain, it can still affect physics at the real energies of the domain from which it originated. 
We illustrate these findings with two simple examples: the scattering problem with the ID oscillating delta-function 
impurity and the static delta- function impurity in quasi- ID wire. 




(29) 



C. Friedel sum rule and Wigner delay time 



V. SUMMARY 
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FIG. 1. (a) Schematic diagram for possible motion of quasibound states in the complex energy plane. E„ and Z)„ denote 

subband threshold energies and domains, respectively, (b) Schematic of the proper Riemann surface R of the domain D„* 
in complex energy plane. The circles and the peaks represent the transmission zeros and the quasibound states (or poles), 
respectively. 
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FIG. 2. Transmission through an oscillating potential for (a) = 0.2, (b) — 0.6, (c) ad = 1.0, and (d) Ud = 1.2 with 
tts = 0, where ad and represent the strength of the oscillating potential and of the static part, respectively (for exact 
definitions see the text). The inset in (a) shows the locations of the poles in complex energy plane for each case. 
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FIG. 3. (a) Trajectory of a quasibound state (the filled circles) as ad increases. The open circles represent the trajectory of 

the quasibound state continued into the neighboring domains on the proper Riemann sheet of Do. The quasibound state of 
the third branch from the top appears from the real axis at ad — 1.869, and collides with sr = 1 and at ~ 2.325 and 
2.935, respectively (see the text for the first and the second branches). Inset: absolute squares of the solution vectors C of 
Eq. (19) at ad = 0.5, 1.5 and 2.5 from the top, respectively, (b) Trajectory of quasibound states as aa increases {as = —0.2). 
At Od = the quasibound state is located at e = 0.8 due to the existence of a bound state, with binding energy eb = 0.2. 
The quasibound state collides with er = at — 0.46. A second branch of quasibound states appears from the real axis at 
ad — 0.525, and collides with er = 1 and at ^ 0.645 and 1.22, respectively. Finally, a third branch appears from the real 
axis at ad — 1.221, and collides with g^j = 1 and at ad ~ 1.475 and 2.035, respectively. 
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FIG. 4. Friedel phases as a function of Floquet energy epi for (a) aa = 0.6, (b) aa = 1.0, and (c) = 0.05 (a^ = 0). (d) 
aa = 0.001 with = -0.2. 
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FIG. 5. The real energies of the trajectory of a quasibound state (crosses) as a function of aa with Os = 0, compared to the 
energies of the local maxima of Wigner delay times (the filled squares) . The insets show Wigner delay times in unit of l/w as 
a function of energy g for aa = 0.4 (left) and aa = 0.8 (right). 
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FIG. 7. Trajectory of a quasibound state for a delta-function impurity in a quasi-lD wire (filled circles) and its continuation 
into the neighboring domain. Inset: real energies of the trajectory of a quasibound state (open circles) as a function of the 
strength of the deltas-function impurity (7 < 0), and the energies of the transmission zeros (crosses). 
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FIG. 9. The real energies of the trajeetory of a quasibound state (open circles) as a function of I7I, compared to the energies 
of the local maxima of Wigner delay times (filled squares). The inset shows Wigner delay time tw in units of 10~^^ sec as a 
function of energy E for 7 = —3.41 feVcm"^, with the local maximum aX E = 3.995. 
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